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The effect of thermal fluctuations on a thin film magnetore- 
sistive element has been calculated. The technique involves 
adding to the basic spin dynamics a general form of interac- 
tion with a thermal bath. For a general anisotropic magnetic 
system the resulting equation can be written as a Langevin 
equation for a harmonic oscillator. Our approach predicts two 
times smaller noise power at low frequencies than the con- 
ventional stochastic Landau-Lifshitz-Gilbert equation. It is 
shown that equivalent results can be obtained by introducing 
a tensor phenomenological damping term to the gyromagnetic 
dynamics driven by a thermal fluctuating field. 



Recently Smith and Arnett have reported that 
thermal magnetization noise may be a serious problem 
for giant magnetoresistive spin-valve heads designed for 
high areal storage density. Using the phenomenolog- 
ical Landau-Lifshitz-Gilbert (LLG) equation with ran- 
dom fields and the fluctuation-dissipation theorem, they 
obtained an estimate for the low-frequency noise inten- 
sity. They considered the OMR head as a single-domain 
anisotropic magnetic particle. 

The goal of this Note is to show that the analysis in 
the framework of conventional stochastic LLG equation 
gives an overstated noise level due to inconsistency of 
isotropic phenomenological damping with basic dynamic 
equations for anisotropic magnetic system. Using normal 
mode approach, we derive a formula for the noise power 
spectrum, which at low frequencies is two times smaller 
than theoretical estimate of Ref . . We also show that 
this result may be obtained in the framework of stochas- 
tic gyromagnetic equation with tensor phenomenological 
damping. 

Thermal noise is a phenomenon that is a result of 
microscopic spin fluctuations due to interactions with a 
thermal bath. It has been customary to include thermal 
fluctuations by simply adding random fields to dynamic 
equations with phenomenological damping, such as the 
LLG equation pj. However the gyromagnetic term de- 
scribes an averaged precession of a large number of spins. 
A direct connection with microscopic physics must be 
found for the phenomenological damping term. 

Lax Q has shown that independent of the specific cou- 
pling to the heat bath, the resulting dynamic equation in 
the vicinity of equilibrium has the form of a damped har- 
monic oscillator driven by a fluctuating force. Here we 



show that for small fluctuations, an anisotropic magnetic 
dynamic system can be transformed into the general form 
analyzed by Lax. We do this by a series of coordinate 
transformations from the original magnetic variables to 
a harmonic oscillator. By back transformation we have 
found a tensor form of the phenomenological damping 
for the magnetization dynamic equation. With this new 
form random fields may now be added to this dynamic 
equation, while properly characterizing the physics of in- 
teraction with a thermal bath. 

Consider thermal fluctuations in a single-domain mag- 
netic particle with the following energy density: 
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Here Vq is the particle volume, and Kj_ are the easy 
(z) and hard (x) axis anisotropy constants, respectively, 
M ={M^,My, M^) is the magnetization (|M| = A/^) and 
is the external magnetic field oriented along the easy 
axis. For a thin film y and z are in the plane. 

For convenience of calculations we introduce the clas- 
sical spin S = — MVb/7i7, where 7 is the gyromagnetic 
ratio. Planck's constant Ti is used as a dimensional pa- 
rameter (it could be another parameter of the same di- 
mensionality, e.g., MsVo/7). All our analysis is pure 
classical and final result does not contain fi. We shall 
describe small oscillations of the magnetization in terms 
of complex variables a* , a, which are classical analogs of 
creation and annihilation operators of a harmonic oscilla- 
tor. This technique is used for a description of spin waves 
(see, e.g., [Q). Complex variables can be introduced by 
an expanded Holstein-Primakoff transformation ^ : 
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The quadratic part of the energy (|l|) can be written as 

E'^'^'^ /n^ Aa*a+{B/'2){aa + aa), (3) 



where A — j{H^ 
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the nondiagonal terms from 
canonical transformation: 



2K^/Ms. To ehminate 
we shall use the linear 



a = uc + vc 
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a = uc' + vc, 
B 

V = — 
\B 



(4) 



A 



1 



where ujq = y/A^ — B"^. Thus, we find the normal mode 
of the system and the energy of an harmonic oscillator 
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^0 = 7V (hH'^ +Ho + K^'>){hI^^ + i/o), (6) 

where loq is the ferromagnetic resonance frequency. 

Using the energy (|^), the dynamic equation for the 
complex amplitude c(i) is: 



dc/dt 
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and is of a form that corresponds exactly to small ampli- 
tude gyromagnetic rotation. Now we consider the inter- 
action of this harmonic oscillator with a reservoir (ther- 
mal bath) 
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which will introduce two related effects: damping and 
fluctuations. Here g^, 5^ are the reservoir variables and 
incorporate specific physical processes. Such interaction 
could have been added to the energy in the a, a* vari- 
ables. However the linear transformation does not 
alter this general form. Following ||], the oscillator 
dynamics for this general form (H) can be represented by 
the following Langevin equation: 



dc/dt = —{iujo + 7])c + f{t), 
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where ry is defined by reservoir parameters and f{t) is the 
random force related to rj. It should be noted that Eq.(||) 
is an exact classical analog of a corresponding quantum 
mechanical equation Note also that phenomenologi- 
cal introduction of dissipation by assigning an imaginary 
part to the resonant frequency coq ± 177 is a common pro- 
cedure in the theory of spin- wave instabilities j7j . 

We can start with Eq.(P) and rewrite it as a stochas- 
tic differential equation with a damping parameter r/ and 
a thermodynamically consistent random field with diffu- 
sion coefficient D ^ {r]/ujo)kBT/h [|, §: 



dc = ~{iujo + v)cdt + VD[dWi{t) + idW2{t)]. (10) 

Here T is the temperature and Wi (t) and W2 {t) are two 
independent Wiener processes. The solution of (10) is 



4^=/ e-^'^°+'^^^'~*''>[dW,{h)+idW2{h)]. (11) 
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Taking into account the rule of Ito's calculus (see, e.g., 
@) dWm,{h)dWnit2) = 5^nHh - ^2)^*1, where A(t) is 
the Kronecker delta function, we can obtain the correla- 
tion function with thermal averaging (...): 

(c*(Oc(0)) = {D/vi) exp(ic^oi - ??*), (12) 
where c*(t) is the complex conjugate of c(i). 



The power spectral density of the complex amplitude 
fiuctuations is defined by 



Sc'c{^) = / (c*(i)c(0))e-^"*dt. 
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Substituting (|T^) into (|T^) and including the correct 
meaning of a damping process for both negative and pos- 
itive time, one obtains 
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Let us consider now an application of our analysis to a 
magnetoresistive head. We shall follow the model [Q of 
a GMR head as a thin film of the volume Vq, uniformly 
magnetized along the z axis of uniaxial anisotropy. The 
transverse anisotropy is defined by shape demagnetiza- 
tion: K± ~ 2ttM^. The power spectral density of ran- 
dom fiuctuations in the voltage dV{t) is defined by 



Svv{(^) = / {SV{t)SV{0))e"^'dt. 
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The voltage deviation for a GMR head is SV{t) = 
Ihias6R(t), which can be expressed in terms of the in- 
plane transverse magnetization ffuctuation 

SV{t)^CvSm,it), Cv = hi.s4r-^T^- (16) 
Thus, Eq.(p^ can be rewritten as 

5yy(w) = CyS,n^jny{l^), (17) 

where 



/ {dmy{t)Smy){0)e''^*dt (18) 
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defines the magnetization noise. 

Since {my{t)) = 0, the deviation Smy{t) — my{t). Tak- 
ing into account Eqs. (p|) and (0), one obtains 



my = -Sy/Sc^-{u + v){c* +c)/V2S. (19) 



Thus 



Smyniy (^) 



{u + vf 
25^ 



{c{t)c*{0))e"^'dt 



{c*{t)c{{)))e''^'dt\. 
Utilizing <^ and (|l|), we have 
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Examination of (^0|) indicates that the damping param- 
eter 77 affects Hnearly the low frequency noise (w <C loq) 
and broadens the resonance peak. However the total 
power is independent of rj. 

In the case of small frequencies uj <^ ujq 
for Hq ^ and At: Ms > Hk- From 



]|) we obtain the 
noise power (as measured for positive frequencies): 



NP ~ 25™„™„(0) ~ rikBT/[TT{jHKMsfVo]. (21) 

The corresponding noise power of 0| (in our notation) is 

NpiSA) ^ 4akBT/[-fiHK)^MsVo], (22) 

where a is the Gilbert damping parameter. To compare 
( pT| ) and (|2^) we need to find a relation between 77 and 
a. This can be done by comparing the dynamic equation 
for the transverse magnetization given in Ref. and the 
corresponding transverse magnetization equation in our 
approach. From [Q we have 

(fmy/dt^ + aAiT^ Msdniy / dt + uJ^my = h[t)^ (23) 

where h{t) is a random "force". Corresponding to (|l 
dynamic equation is: 



d^niyl dt^ + 2rjdmy/dt + ijJ^niy = h{t). 



(24) 



From (|23|) and (24), we obtain the following relation: 
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Substituting (|25|) into ( |22|) , one can see that at low 
frequencies our theory (pTl) predicts 2 times smaller 
(3 dB smaller) noise power than the result of [Q: 
NP/NP'^^^^ = 1/2. It should be also noted that the con- 
ventional approach with stochastic LLG equation gives 



different noise frequency dependence than our result (20) 



Now we demonstrate the inconsistency of the con- 
ventional phcnomenological damping term with basic 
equations for a harmonic oscillator and show how this 
isotropic damping can be replaced by a tensor damping 
term to obtain the above results. 

From the LLG equation 

dM a"/ 

— = -7M X Heff - X (M X Heff ) (26) 

one can derive the following linearized equations: 
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On the other hand, we can obtain the linearized equation 
for TOa, and my from the normal mode equations (Eq. (ITT 
without the thermal term) 



{d/dt + ri)c = —iujqc, {d/dt + ifjc* = iujqC* 



(28) 



Substituting c = ua ■ 
taking into account ( 



wa*, c* = ua* — va into ( pSD and 
]) and m — ~S/S, we obtain 
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By comparison with ( P7| ) we see that the nondiagonal 
terms coincide. The diagonal terms in ( ^ ) and ([2S|), 
responsible for relaxation, differ from each other. 

In order to have the magnetization equations (p9[), one 
should start with the magnetization equation containing 
a tensor damping of the form: 



^ = -7M X Heff - 7^ X [a .(M X Hes)], 
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For small oscillations about equilibrium, the a^z compo- 
nent does not enter into the equations. The conventional 
LLG equation can be used only for a uniaxial system, 
where H^'' = and the tensor is effectively a scalar. 
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